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ABSTRACT 

I use a 1-D Lagrangian code which follows both a gaseous and a dark component to 
study the effects of a metagalactic UV background on the formation of dwarf galaxies. 
In the case of an ionising background consistent with the existing observational con- 
straints no effect has been found in halos with circular velocity above 60 km s _1 . The 
effect is very mild in halos with 40 < V c < 60 km s -1 , and the reduction in the total 
gas mass which cools below 8000 K never surpasses 30 %. Below the circular velocity 
40 km s^ 1 there is a sharp transition towards the case in which the gas does not cool 
but it is gravitationally confined by the potential well of the dark halo. Halos with a 
virial temperature below the equilibrium photoionisation temperature cannot confine 
the gas. 
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In the standard picture of structure formation smaller ob- 
jects collapse first, and merge later on to form bigger sys- 
tems. As the gas falls into the potential well of the form- 
ing halos it passes through a shock in which its bulk ki- 
netic energy is converted int o thermal energy. In Paper I 
(Forcada-Miro & White 1995) we studied how the hot pri- 



mordial plasma can lose its thermal energy via radiative 
cooling. The gas cools more efficiently at higher redshifts, 
due to the bigger densities, and in smaller halos, due to line 
cooling. The result is then too many dwarf galaxies, which 
causes the faint end of the luminosity function to be too 
steep compared to observations (c.f. White & Rees 1978; 
Kauffman, White & Guiderdoni 1993; Cole et al. 1994). 

Several possible solutions to the overcooling problem 
have been contemplated in the literature. A first possibil- 
ity is energy feedback from the stars which form out of the 



cold de nse gas sitting at the centre of the halo. Dekel fe Silk 



(1986) argue that supernova explosions could produce a big 
gas loss in halos up to a circular velocity V c ~ 100 km s _1 . 
However, the physics of star formation and how stars in- 
teract with the surrounding interstellar medium is still very 
poorly understood. Another possibility is a higher merging 
rate, but then one could run into the problem of thicken- 
ing too much the discs of spiral galaxies (Toth & Ostriker 
1992; Navarro et al. 1994; Velazquez & White 1998). A third 
possibility was first suggested by Ikeuchi (1986) and Rees 
(1986) as a model for the Ly-a clouds: a UV background 
could prevent the primordial gas from cooling and sinking 



towards the centre of small halos. Efstathiou (1992) pointed 
out that an ionising background compatible with the ex- 
isting observational constraints could reduce the predicted 
number of dwarf galaxies. 

These three possibilities are not exclusive, and the final 
answer to the overcooling problem is probably a combina- 
tion of all of them. In the present paper I explore the role 
of the UV background in the process of galaxy formation. 
The strongest piece of evidence for the presence of a meta- 
galactic ionising backgr ound is the absence of the so-called 
Gunn- Peterson trough (Gunn & Peterson 1965) in the spec- 
tra of quasars (Steidel , 
Bahcall et al. 1991; 



Sargent 1987; Schneider et al. 1989; 
Webb et al. 1992). Such a low Ly-a op- 
tical depth could also be explained if there was no inter- 
galactic medium (IGM), but this is quite improbable since 
the process of galaxy formation is not thought to be 100 % 
efficient. In addition, the detection of intergalactic ionised 



helium absorpt ion by the Hubble Space Telescope ( lackob- 



sen et al. 1994) has not only confirmed the existence of a 



highly ionised IGM, but also the presence of a substantial 
amount of helium as predicted by the Big Bang nucleosyn- 
thesis. 

The present paper is organised as follows: in Section 2 
I briefly describe the model assumed for the ionising spec- 
trum; in Section 3 I present the numerical simulations done 
with the spherical Lagrangian code described in Forcada- 
Miro & White (1998); the results obtained are discussed in 
Section 4. In all cases the Hubble constant has been taken 
to be H — 50 km s _1 Mpc -1 , and the baryonic content of 
the universe is Ob = 0.05. 
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2 MODEL FOR THE UV BACKGROUND 

In the simulations presented in this paper I have assumed 
an ionising spectrum which is a simple power law, 



In the numerical simulations presented in the following 
section we refer as a fiducial model to the case with a UV 
background which has J-21 = 1, a = 1.5, and Zi on = 5. 



J„ = J Vh {z)[ — 



(1) 



where a is the spectral index, and J„ L is the mean intensity 
of the radiation field at the Lyman limit. The ionising back- 
ground is switched on at a redshift Zj on , and its subsequent 
dependence on the redshift is taken as in Efstathiou (1992), 



j„ l 0) = 







l+[5/(l+z)] 4 



Z Z> Z[ on 
Z Zion 



(2) 



where J-21 is in units of 10 21 ergs cm 2 s 1 



Hz 



The proximity effect provides an estimate of the inten- 
sity of the UV background at the Lyman limit ( Lu et al 
199ll), 



log J„ L = -21 ± 0.5 for 1.7 < z < 3.; 



(3) 



where J„ L is in units of ergs cm -2 s^ 1 Hz -1 sr _1 . More 
recent determinations based on high resolution spectra do 
not find also significant evidence for a change in J^ L in 
the redshift range 2 < z < 4, but they seem to favour 
a smaller value ~ 3 x 10 -22 ergs cm -2 s _1 Hz -1 sr _1 
(Willinger et al. 1994; Fernandez-Soto et al. 1995; Gial- 
longo et al. 1996; Cooke et al. 1997). At lower redshifts 
Jv L is still poorly constrained. Kulkarni & Fall (1993) 
have analysed HST observations, and obtain a value ~ 



All these estimates 



6 x 10~ 24 ergs cm" 2 s" 1 Hz" 
are determined assuming a power-law distribution in red- 
shift of the Ly-a clouds, and it could be the case that such a 
distribution is not adequate. In the simulations I have con- 
sidered a stronger background, with J-21 = 1, and a weaker 
background, with J-21 = 0.1. The existing observational 
constraints fall in between these two cases. 

A possible source for the ionising radiation is quasars, 
which have a spectrum that in the extreme ultr aviolet is 
well app roximated by a simple power law a « 1.8 (Zheng ei 



al. 1997). However, the spectrum of quasars is too hard to 
be in agreement with Hell observations (Madau & Meiksin 
1994; Savaglio et al. 1997). Also if the population of quasars 
has a strong decline at z > 3 the emitted UV radiation 
i s insufficient to ionise the medi um to the required level 
( [Miralda-Escude fc Ostrikcr 1990| ). Recent redshift surveys 
seem to indicate that the decline in the quasar population 
is not so strong (Boyle 1991; Irwin et al. 1991), and then, 
t he ionising flux provided by the q uasars might be sufficient 
(Miralda-Escude & Ostriker 1992). Another possible source 



of UV radiation is young stars in galaxies, which only emit 
below the Hell ionisation edge. Therefore models including 
the contribution from both quasars and starbursting galax- 
i es provide a softer spectrum than models with quasars only 
(Miralda-Escude & Ostriker 1990). Also the first generation 



of baryonic objects in which stars could form might emerge 
at redshifts z gj 10 — 50 thanks to the formation of molecu- 



lar hydrogen ( Haiman fc Loeb 1997 ). Therefore, I have con- 
sidered both a hard background, with a = 1.5, and a soft 
background, with a = 5, and the ionising flux is switched 
on at two different epochs, Z[ on — 5 and Zi on = 10. 



3 NUMERICAL RESULTS 

3.1 Effects on halos with circular velocity 

V c > 40 km s" 1 

In the presence of the fiducial UV background the be- 
haviour of the radiative shocks differs significantly from the 
behaviour discussed in Paper I up to a circular velocity 
V c ~ 60 km s _1 . Shortly after the ionizing background is 
switched on at Zi on = 5 there is a short period of time in 
which the shock wave moves outwards very fast, approaching 
the nonradiative behaviour (Figure 1). This effect is stronger 
as one goes towards lower circular velocities. Notice also that 
the oscillations with growing period and amplitude which 
appeared in the simulations with no ionizing background 
are not present in these calculations. 

In Figure 2 I compare the effect of varying some of 
the parameters of the model in the case of a halo with 
circular velocity V c — 50 km s _1 . If the background is 
switched on at a higher redshift (i.e. z- lon = 10) the shock 
approaches the non-radiative behaviour at an earlier time, 
but the subsequent evolution is very similar to the fiducial 
model. If the intensity of the ionizing radiation is decreased 
(i.e. J-21 = 0.1) some small oscillations appear in the evo- 
lution of the shock wave. The most drastic effect is obtained 
by softening the spectrum (i.e. a = 5), since the oscillations 
with growing period and amplitude are present again. 

On the other hand, the evolution of the total mass which 
has been able to cool below 8000 K does not differ so much 
from the evolution when J-21 = 0. For redshifts 3 < z < 5 
the presence of the UV background does not affect at all 
the ability of the gas to cool inmediately after being shock- 
heated. It is not until redshifts z < 3 than the mass of cold 
gas is smaller in the presence of the fiducial UV background, 
but in any case, the reduction does not surpass 30 % (Fig- 
ure 3). Ionizing the gas at a higher redshift or reducing the 
intensity of the background does not affect practically the 
amount of cold gas. However, in the case of a softer spec- 
trum, the evolution of the cold mass is closer to the evolution 
in the absence of the ionizing background. 



3.2 Effects on halos with circular velocity 

V c < 40 km s- 1 

I have simulated as well a halo with circular velocity V c = 
30 km s _1 . In this case the evolution of the system is more 
complicated and it is better understood by looking at the 
radial profiles of the velocity, density and temperature at 
different redshifts. In Figure 4 one can see that an outflow 
forms, which is driven by the photoionizing radiation. As the 
system evolves the photoionization-driven wind gradually 
dies out and the cloud starts growing via cosmological ac- 
cretion again. The result is then a cloud in quasi-hydrostatic 
equilibrium in which hot gas is confined by the gravity of the 
dark minihalo. Shortly after switching on the ionising back- 
ground the gas cannot cool and settle towards the centre of 
the system anymore. 
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As in the case for bigger halos, these results are not very 
sensitive to reionizing the gas at an earlier time. If the flux 
of the UV background is reduced the wind is much weaker 
and the gas has some residual infall. In the case of a much 
softer spectrum the ionizing radiation cannot drive a wind, 
and the gas is able to fall in the potential well of the halo 
and cool efficiently. 



estimate the temperature of the intergalactic medium by 
nearly an order of magnitude in the early stages of reion- 
isation. This is because the photoionisation time is much 
shorter than the recombination time just after the back- 
ground is switched on, and the ionising radiation can de- 
posit energy more efficiently in the plasma. At later times, 
radiative recombination becomes important, and both cal- 
culations converge. 



4 DISCUSSION 

There are two different ways in which the UV background 
can influence the ability of the primordial plasma to cool and 
collapse. One of them is by increasing the ionised fraction 
of the gas, and thus, reducing line cooling. The other is by 
heating the IGM by producing electrons with excess kinetic 
energy. 

The efficiency of the ionising radiation to modify the 
emissivity of the plasma depends on the hardness of the 
spectrum and the density of the plasma (Figure 5). If the 
spectrum is hard both hydrogen and helium are strongly 
ionised, and the cooling function does not have such a neg- 
ative power-law dependence on the temperature, but a soft 
spectrum does not affect the ionisation state of helium. This 
explains why in the case of a hard spectrum the shock wave 
is not oscillatory unstable, but the oscillations are present if 
the spectrum is soft. On the other hand, collisional ionisa- 
tion rates are proportional to p g 2 , since they are two-body 
processes, while photoionisation rates are only proportional 
to p g . In this way, if the density is low enough, line cooling 
is heavily suppressed, and the emissivity of the plasma is 
reduced up to temperatures T ~ 10 6 K. However, at the 
post-shock density, which is o n average of the order of 12 
times the background density ( Forcada-Miro fc White 1998 ), 
the emissivity is only reduced up to a temperature few times 
10 s K. This explains why no significant effect is obtained for 
halos with a circular velocity V c > 60 km s _1 . The gas is also 
able to cool efficiently enough in halos with circular velocity 
down to V c ~ 40 km s~\ and thus, there is not so much 
difference in the total cold mass. However, the reduction in 
the emissivity of the plasma increases the cooling length and 
the shock wave moves outwards faster until the free-fall time 
is of the order of the cooling time (see Paper I). 

There is a temperature T eq , known as the photoionisa- 
tion equilibrium temperature, at which the energy losses A 
equal the heating due to the photoionising radiation V. As 
one goes towards smaller halos the virial temperature is only 
slightly bigger than the equilibrium temperature, which in 
the fiducial case is T cq ~ 3 x 10 4 K. Consequently, in ha- 
los with a circular velocity in the range 25-35 km s _1 the 
gas is still gravitationally confined, but the heating due to 
the photoionising keeps it at a temperature of the order of 
T eq . Such systems were proposed by Ikeuchi (1986) and Rees 
(1986) as possible candidates for the Ly-a forest. This ve- 
locity range is in agreement with Meiksin (1994), and is very 
much independent of Zi on and J-21. However, in the case of a 
softer spectrum the equilibrium temperature is only 10 4 K, 
and the Ly-a clouds would have a characteristic circular ve- 
locity V c ~ 20 km s -1 . 

Finally, I would like to point out that a time-dependent 
calculation of the ionisation state of the plasma is crucial, 
since calculations assuming equilibrium abundances under- 



5 CONCLUSIONS 

The qualitative conclusion I draw out of the simulations is 
that the metagalactic UV background has a very mild ef- 
fect on the formation of dwarf galaxies. The way the ionis- 
ing radiation affects the collapse of small objects is more 
by heating the pre-shock plasma than by modifying the 
cooling function. Our results are only marginally in agree- 
ment with the spherical collapse calculations of Thoul & 
Weinberg (1996), which indicated a reduction ~ 50% of 
the cooled mass fraction in halos with V c ~ 50 km s , 
and a negligible effect above V c ~ 75 km s . They find a 
stronger effect probably because they use a slightly harder 
spectrum (i.e. a = 1), and as I have showed the effect of 
the UV background is very sensitive to the spectral index. 
On the other hand, the three-dimensional simulations of 
Navarro & Steinmetz (1996) show a mild effect of the meta- 
galactic background all the way up to normal galaxies (i.e. 
Vc ~ 200 km s _1 ). As the authors themselves suggest, this 
must be due to the fact that a fraction of the gas which 
falls onto the potential well of a bigger halo collapsed at an 
earlier time onto smaller halos. 

In any case, all the theoretical work seems to confirm 
that a UV metagalactic background is far from being enough 
to solve the overcooling problem, since the luminosity func- 
tion semi-analytic models predict differs from observations 
at circular velocities V c < 100 km s _1 . 
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Figure 1: Evolution in time of the shock wave position (solid line) and the cooling wave position (crosses) for halos with 
40 < V c < 60 km s" 1 in the presence of the fiducial UV background. Long-dashed lines correspond to the adiabatic evolution of the 

shock. 
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Figure 4: Radial profiles at different redshifts for a halo with V c = 30 km s _1 in the presence of the fiducial UV. Top: Velocity. 
Middle: Density. Bottom: Temperature. The variables have been scaled as in Paper I. 
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Figure 5: Cooling and heating functions for a primordial gas. In all cases Yp = 0.24 and J— 21 = 1. Compton cooling is not considered. 
Left panel: Effect of density. Thin lines correspond to p g = 5.864 X 10 — 29 gr cm -3 . Thick lines correspond to a density 200 times 
bigger. In both cases the spectral index a = 1.5. Right panel: Effect of spectral index. Thin lines correspond to a = 1.5. Thick lines 

correspond to a = 5. In both cases p g = 5.864 X 10 — 29 gr cm -3 . 
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